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Preamble 


Retrospective recollection and inspirational outlook 


I have been asked by our publishers to contribute a preamble to the Cumulative Index of 
Volumes 1-100 of our Journal. It is indeed hard to believe that the first number appeared as 
recently as June 1972, just short of twenty years ago. More astonishing is that in the first two 
years 1972 and 1973, only one volume could be published in each year. Switching, however, 
forward to the year 1993, I can inform our readers that 9 volumes will appear. This 
demonstrates the stormy development that has taken place in the evolution of the main 
thematic topic of our Journal: the exploration of the computer simulation of complex physical 
and technical problems. We should remember, however, that the revolution in our thinking 
processes and their application to problems of industry was initiated as far back as in the 
closing years of World War II under the impact of the first appearance of electromechanical 
computing devices in the United Kingdom and the United States and the conception of 
alternative new formulations of mechanics. But it is equally important to recollect that 
practically all efforts in the first twenty years or so in engineering computing science (say until 
1963 or so) were directed towards the solution of problems in solid mechanics, e.g. those 
arising in plates, shells, wings, fuselages, etc. But the real revolution was actuated — and this is 
indeed in general unknown — through the then extremely severe problem of the structural 
safety of swept-back wings which, for aerodynamic reasons, had to be provided for fighter 
aircraft under jet propulsion in the last two years of World War II; these combat aircraft were 
designed for high subsonic speeds. Initial thoughts on a suitable theoretical framework, far 
removed from the clumsy and inefficient relaxation methods of the late thirties, were first 
generated in 1944 and led to the birth of the so-called matrix displacement method and the 
conception of a triangular membrane element (a true deus ex machina for swept-back wings, 
now a standard component of the finite element method) as well as the first matrix software 
codes in the UK and shortly afterwards (and no doubt independently so) in the USA. 
Germany did not then participate in this mental revolution. This preoccupation with structural 
problems, novel and significant as they were in aeronautics, was dictated by the slow speed 
and restricted memory of the first digital computers. On the other hand, major breakthroughs 
in the non-linear analysis of structures were already achieved in the late fifties by the inventive 
definition of a geometrical stiffness — which proved highly successful in instability and large 
deflection problems — and the concept of initial load or initial strain, which made possible the 
analysis of structures possessing non-linear material properties, e.g. plasticity, viscosity. 
Effective solutions to such problems were developed in the late fifties in the UK and the USA. 
An excellent presentation of the state of the art was given in 1965 at the exhilarating 
Dayton-Ohio Conference on Matrix Methods in Structural Mechanics organised by John S. 


Elsevier Science Publishers B.V. 
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Przemieniecki, which proved a landmark in the development of computer mechanics in the 
solid domain. 

But already in the mid-sixties, the necessity of developing CFD beyond the then available 
finite difference schemes became increasingly urgent and led to the development of techniques 
in the subsonic regime mainly based on FEM and BEM. This evolution was activated through 
the appearance of faster and more powerful digital computers. In particular, machines like the 
CD 6600 and CD 7600 greatly accelerated the developments and generated a stream of 
publications towards our Journal. 

In the preceding paragraphs, our thoughts were mainly focused on the past history of the 
Journal. However, a scientific Journal lives and prospers only through daring excursions into 
novel research activities. In fact, our Journal published in the past years a number of 
distinguished papers on reentry aerodynamics of space shuttles — but more of this later. Other 
aerodynamic problems were also discussed. 

The preceding remarks were intended to demonstrate how novel fields of research 
necessitate a continuing realignment of editorial policy. Thus, it may be pertinent to point out 
that the FEM has conditioned a mental revolution in the transformation of engineering 
sciences but may at the same time have reached the limits of its applicability in certain areas 
like CFD, in particular in the presence of turbulence. New simulation techniques and concepts 
are urgently called for. This is a major challenge to authors and editors alike. 

Indeed, we are happy to report that Computational Fluid Dynamics has been a rapidly 
growing field of scientific and engineering interest over the last twenty years or so. We 
mentioned above that such developments were closely associated with the availability of even 
faster and more powerful digital computers. In this evolution, steady-state methods based on 
potential or stream functions were evolved, using first Cartesian and subsequently body fitted 
grids. Such formulations were found useful in the application of the finite element method to 
the elliptic, subsonic regime. Moreover, finite difference and finite volume techniques were 
developed in the unsteady regime but were essentially restricted to a one-dimensional regime. 
As a result of the continuously increasing computer power, extensions of FEM permitted the 
steady and unsteady solution of two-dimensional problems using initially the Euler equations 
and subsequently also the Navier-Stokes equations. The finite element method has in 
particular been applied to incompressible and compressible viscous flow. Modifications of the 
Navier-Stokes equations made it possible to include effects of turbulences via a system of 
time-averaged equations. 

The revolutionary appearance of vector computers in the eighties at last gave us the 
computing tools to solve two- and three-dimensional problems using increasingly fine grids. 
Initially, these architectures favoured explicit time-stepping schemes which proved easy to 
vectorise. Subsequently, new strategies for the vectorisation of solution techniques also 
improved greatly. 

Although first attempts for a direct numerical solution of turbulent flows have been 
initiated, even the most efficient present day computers are not in a position to compute a 
fully turbulent flow within a prescribed simple geometry. New simulation techniques and the 
evolution of parallel computer architectures possessing sophisticated software are required to 
find a solution to such complex problems. Although the FEM achieved triumphant results for 
laminar regimes and facilitated the analysis of reentry flights by the incorporation of real gas 
effects, its scientific limits are clearly discernible as mentioned above. Of course, we have to 
recognise again that the FEM has achieved very fine results for compressible flow problems in 
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hypersonic aerodynamics due to its inherent capabilities to handle unstructured meshes and 
adaptive mesh refinements. Several techniques have been evolved and applied with signal 
success in the past years. From the mathematical point of view, the FEM undeniably provides 
a secure philosophical foundation, but in the presence of compressibility, it does not yet 
permit a consistent application of this fundamental advantageous framework. 

Focusing our attention once more on the hypersonic regime, we have to point out another 
conceptual limitation of present day CFD techniques. Thus, we note that the current FEM 
applies in the hypersonic regime a linear approximation within each element and consequently 
provides a solution which does not differ essentially from that of the finite volume technique. 
At the same time, finite volume and finite difference methods compensate their deficiency in 
the realm of unstructured meshes and adaptive mesh refinement as well as in the simplicity of 
boundary conditions by offering in contrast a high sophistication through so-called high 
resolution schemes. Although successful attempts have been made, the consistent application 
of such strategies to FEM is by no means straightforward. 

To appreciate the complexity and magnitude of the simulation process in the extreme case 
of hypersonic aerodynamics, consider, for example, a reentry of a space shuttle, including, 
moreover, dissociation and reassociation processes of the air. Thus, at a first rough estimate at 
M = 25 at an altitude of 75 km, a typical grid of a finite element scheme possesses 589 865 
unknowns (symmetrical re-entry manoeuvre). This number is, however, increased to a 
multiple if the dissociation process as well as the viscous transport processes in the air are 
included. Let us also admit that this computing simulation is still incomplete if we do not also 
include the effects of turbulence at least over part of the air stream. Moreover, the above 
considerations apply strictly to the case of stationary aerodynamics whilst we now have to 
consider the necessity of a numerical analysis of instationary aerodynamics in which remark- 
able progress has already been achieved. It is clear from this account that such major tasks 
cannot be efficiently solved on existing vector supercomputers like the Cray 2 and the NEC 
machines and their immediate successors but demand parallel computational configurations of 
the highest sophistication as initiated in the USA, Japan and Europe. It must be admitted, 
however, that the efficient and philosophically sound design of parallel software is still in its 
infancy. 

Two further examples serve to underline the necessity of such developments. The first 
refers to superplastic forming processes. The aerospace industry is striving increasingly to 
produce high-quality component parts, such as satellite tanks, using superplastic forming. On 
commission from Aerospatiale of France, a novel FEM methodology was evolved in the early 
eighties specifically for the simulation of superplastic forming processes and has been applied 
successfully to a number of aerospace structural components. More recent developments 
realise superplastic forming processes by means of numerical simulation and significantly so by 
the optimisation of the forming pressure to achieve a minimum processing time as well as the 
determination of the initial thickness profile of a given sheet for a net-shape forming process. 
Such simulation processes and other related investigations like those of thermally coupled 
forging processes may require 50 hours of CPU time on a Cray 2, the turn around time being 
four or five times larger. The computational effort is similar in the aerothermodynamics of a 
complete — relatively simple — flying configuration like the Hermes space shuttle. And yet the 
numerical simulation of combustion processes requires even higher computer capacities, which 
can only be achieved through highly parallelised novel computer configurations as they are 
being developed presently. 
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The aforementioned combustion processes possess an increasingly significant and critical 
environmental significance. In fact, we are expected to attain higher combustion temperatures 
in order to improve the combustion processes. The aim is not only an increased energy yield 
but also a drastic reduction of pollutant emissions. The Swiss firm ELCO Energy Systems, 
which, together with its European subsidiaries, covers a significant part of the heating and 
industrial combustion market in Western Europe, sees an obligation to tackle this problem in 
a fundamental and long-term way. Moreover, it has been recognised that only novei physical 
and numerical techniques can lead to major innovative results. The set task includes the 
realistic simulation of combustion processes in turbulent three-dimensional gas and two-phase 
flows and the optimum conception of the combustion chamber geometry as well as the 
thermally highly loaded materials. A novel simulation programme is being designed specifical- 
ly for the exploitation of compact high-performance parallel computers. 

The above topic brings us to the evolution of a future computer-based materials science 
and the provision of new materials like ceramic coatings. Such design tasks can again only be 
realised with the aid of novel computer simulations. Increased performance of turbine aero 
engines, for example, can only be attained by higher gas temperatures which exceed the 
thermal stress capacity of the current turbine blades. To protect such blades, ceramic coatings 
have to be sprayed onto the surfaces by means of a plasma beam. Due to the production 
process, the result is a complex material layer structure with porosities, the thermomechanical 
resistance of which is unknown; yet it is decisive to elucidate and quantify its nature in order 
to safeguard the service of the component. The same problem also arises in the automotive 
industry where ceramic-coated lightweight pistons are to be exposed to higher temperatures in 
the combustion engines. Especially the computer simulation of the coating process on the one 
hand and the determination of the thermomechanical behaviour of the coating on the other 
hand necessitate the development of novel computational methods which do not rely only on 
FEM, but also on alternative simulation (fractal) techniques. Such methods have to account 
for the statistically random distribution and growth of flaws and fissures as well as the 
evolution of residual stresses. 

Our preamble can only spotlight a few problem areas which have to be solved effectively 
and economically through the evolution of new simulation techniques and modern parallel 
computer devices if supported by highly imaginative software. To reflect and propagate such 
progress remains the noble task of this Journal. We can only hope that a new anniversary like 
a forthcoming 200th volume will demonstrate that, within a few years, some of our 
expectations have been fulfilled. 

We can only hope that our condensed account will encourage old and new authors alike to 
partake in this eternal race to a momentary perfection. Inevitably, our survey reflects to a 
certain extent personal preferences and experiences of the editorial board. In this sense, we 
selected CFD and pointed out some of its intriguing major problems which demand and have 
not found yet an advanced application on present day and future computing devices. Similar 
considerations apply to the computer evolution of materials science. These two domains of 
CFD and material science associated with an inventive application of parallel devices will 
certainly occupy this Journal over the next 10 years or so. These scientific aspects also 
demonstrate how the interests and unsolved problems which occupy the mind of authors and 
editorial board alike change (and have to change continuously) the scientific orientation of a 
journal as long as the editorial board is determined to remain at the top of the league. 
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However, there is still another major field which will occupy us: this concerns the 
explosively developing theory of chaos and its scientific and industrial applications. That this 
necessitates a highly sophisticated computer simulation we need not even emphasize. The 
theory of chaos has emerged with such explosive force and imaginative vigour following the 
pioneering findings of E.M. Lorenz who wrote a classic paper in 1963 on the evolutions in a 
Rayleigh—Bénard layer selected to simulate the meteorological response of the atmosphere. 
This publication demonstrated first that deterministic systems may respond chaotically, this 
being expressed by the appearance in the phase space of the now famous strange attractors. 
Indeed, this paper initiated a new era. In fact, many leading physicists but also engineers are 
of the opinion that the theory of chaos is the science of the next century. Clearly, a journal 
like ours has to direct its attention to this and other novel manifestations of computer science. 
There is also no doubt that the understanding of chaos will be one of the most significant keys 
to the unravelling of the mysteries of turbulence. Such research inquiries will lead to new 
computer simulations, for example in CFD, including ab initio the effects of possible 
turbulence. We hope that our authors will contribute energetically to such novel scientific 
findings and publications. 

There is no doubt that the FEM proved superior to other techniques in adapting to 
complex configurations and boundary conditions in physical problems; at the same time, 
however, the FEM technique involves higher computing times than, for example, the finite 
difference technique. Thus, we have to realise the verdict of history “Ta mavta pet” and the 
consequent necessity of reassessing past attitudes. We wish our readers the greatest success in 
helping us to unravel some of the unsolved problems of engineering science. 

We on the editorial board have to accept that this necessitates very imaginative efforts, not 
only on the scientific side. We have to admit that without the devoted support of the members 
of the confraternity of our publishing editors and of the devoted work of the desk editors of 
our Publishing House in Amsterdam, such success and speed of action cannot be maintained 
over long periods of time. 

In the course of the 20 years or so, we had the support of an array of superb publishing 
editors like the then director Wim Wimmers of North Holland and his successors Einar 
Frederiksson, Arjen Sevenster, Arie Jongejan, Bas van der Hoek, Hans van der Nat and 
Mark Eligh. 

To all those, we express our deeply felt gratitude for their devoted and imaginative support 
and for having guided the editorial board with such steady hand. Particular thanks are 
expressed to the magnificent Ineke Kolen who with unflappability and a penetrating mind 
supports us smoothly in our difficult task. Wim Wimmers, who initiated us in the art of journal 
publishing, demonstrated to us how such a venture can lead to a steady success in science and 
engineering and also be economically viable. A special expression of gratitude is also extended 
to the last two members of the publishing editors Hans van der Nat and Mark Eligh who 
understood our problems, helped in many difficult situations and also became good friends of 
the editorial board. 

Ultimately, we remember our guardian angel, the editorial secretary Marlies Parsons, who 
guided us with Swiss clockwork precision through many stormy days. 


John Argyris 
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